Coulomb Screening of 2D Massive Dirac Fermions 
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A model of 2D massive Dirac fermions interacting with a instantaneous 1 jr Coulomb interaction, 
is presented to mimic the physics of gapped graphene. The static polarization function is calculated 
explicitly to analyze screening effect at the finite temperature and density. Results are compared 
with the massless case . We also show that various other works can be reproduced within our model 
in a straightforward and unified manner. 
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Over the past several years, the physics of graphene has attracted considerable interest, both theoretically and 

experimentally [1]. Graphene, due to its two dimensional hexagonal lattice structure, has unique linear energy spec- 
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trum near the Dirac points of Brillouin zone. Because of the unusual energy band dispersion, many electronic 



q 



properties in graphene exhibit significantly different behavior from the conventional 2D systems, for example, half- 
integer QHE(Quantum Hall Effect) Q|. 

In graphene although the motion of electrons are fixed on the 2D plane, their interactions still show 3D Coulomb's 
i-^j , law, for the electric field lines cannot be confined on 2D. While most of the early work were based on massless Dirac 
q ■ fermion model, recent work has shown that the massive case can also be created[2j. 

i For massive fermions, it is equivalent to the opening of a gap in the electronic spectrum in condensed matter 

| physics. So the similarity between graphene and QED2+1 is obvious, and many results were obtained by explor- 
ing this correspondence[3j. While the main difference is that Graphene lacks of Lorenz invariance due to its nearly 
■ instantaneous coulomb interaction. 



fT} ' Besides, it should be noted there has been many papers 
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on the Coulomb interaction in gapless graphene. 



Ref.PlQl gave the polarization function for zero temperature gapless graphene at finite density. Ref.[(| dealt with 

2D Coulomb- interacting massless Dirac fermions and calculated the specific heat at finite temperature, while Ref.[7] 

• • 1 generalized it to the finite density case. As for the massive Dirac fermions, there are also many work on it. However, 
^ ■ 

most of them (gapless or gapped graphene) were highly succinct, used different models and did not give the calculations 
?H ' of polarization function in detail. 



In this paper, we consider a model of two-component Dirac fermions interacting through a three dimensional 
instantaneous Coulomb interaction, and calculate the polarization functions at finite temperature and finite density 
using finite temperature field techniques [8]. In fact quantum field theory at finite temperature or density is usually 
applied to study cosmology and astrophysics But its 2D spatial case has not been observed in nature. So the graphene 
provides a wonderful platform for establishing QED2+\. We also show various works on this topic can be connected 
within our model in a straightforward and unified manner. 
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II. MODEL 

Our starting point is a model of 2 + 1 dimensional two-component Dirac fermions mediated by a three dimensional 
Coulomb interaction at temperature T. The action of the system S is given by (h = 1) 

f 13 f 1 f N - 

S(i/j,ip,(p) = / dT{ d 3 x — \d l ip(x,T)\ 2 + d 2 x^2'ip s (x,T)[dr+vcr-p + m a 3 
J o J n J s=1 

-(j, + ie(p(x,T)]i> a (x,T)} (1) 

Here, j3 — 1/T and \i is the chemical potential. The fields ip s are two-component fermion fields, subscript index 
s stands for different species of fermions with N = 4 due to spin, valley degeneracies in graphene. The vector 
a = (oi, &i) and Oi, i = 1, 2, 3 are pauli matrices; oo = 7 is the 2x2 identity matrix(omitted for simplicity), (p is the 
field that mediates the Coulomb interaction, nothing but the time component of the electromagnetic field. The 3D 
space integral of the action describes the kinetic term of the scalar field, while the remaining terms describe kinetic 
term for the fermion fields and their interaction with the scalar field. In addition, we put v = 1 bellow and only 
restore it if necessary. 

Correspondingly, the Green's function for the free Dirac fermions 

fc + <t ■ k + ma 3 

G o(fc) = — - — -o — (2) 
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Using more symmetric three-momentum notation, (qo , q) = (icJi, q), (fco, k) = {iu) n + /i,k),fc = k 
(2n + 1)tt//3,iji = 21?:/ (3 etc., then the polarization function in the random-phase approximation (RPA) 



is 



f rf 2 k fc (fc + g0 ) + k(k + q) + m 2 

where u(q) = 27re 2 /«|q| is the 2D Fourier transform of the 3D Coulomb interaction. Following similar consideration 
in Ref.[10] we can divide II into two contributions, the vacuum part and matter part. 



II — Hvac + ^matter (4) 



lim IT = U vac (5) 

The vacuum polarization function, both massless and massive case, can be calculated by dimensional regularization 
approach as in Ref . Q ■ Here, we can reproduce their results within our model 



Iql 2 , m 1 ,„ 4m 2 . , q . 

U vac (q ) = -^- {^ + - 1-— arctan (JL } 6 
^ 2g g^ 2m 
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III. POLARIZATION FUNCTION FOR MASSIVE DIRAC FERMIONS 



The sum of fermion Matsubara frequencies can be performed in a standard manner 



1 ^ fc (fc + g ) + k • (k + q) + m 2 



0. 



With 



{K-El)[{k + q l)-El +(l ] 

= -—. <t dzh{z)g{z) 
2m f 

= / dzhiz)- tanhCfe /2) / dzh(z)-t&nh(f3z/2) 

2tt« Jioo-e 2 2m J-ioo+e 2 

i p— ioo — e ii i pioo+e i i 



, _ (z + n)(z + (1 + g ) + k ■ (k + q) + m 2 



[(z + tf-EfiKz + v + qft-Ei 



k+qJ 



g(z) = | tanh(/3z/2) 



£k = \/k 2 + to 2 

For function g(z) has simple poles at z = zu;„, the sum emerge as the integration of the product /ig along a suitable 
path in the complex plane. We can divide the above expressions into two parts ir vac , ir ma tter, which describe the 
vacuum and matter's contributions respectively. 



W = ^zr / dz[h(z) + h{-z))/2 (8) 

-ioo 
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Where we have defined 



Y pioo+e ^ ^ p — ioo — e ^ 

^ = -^iL 00 J zh{z) ^Ti + ^iL-, dzh{z) ^r-i 

1 £ k (£ k + g )+k(k + q)+m 2 1 
2£ k (£ k + go ) 2 - £k+ q e^+z*) + 1 

1 g k (£ k -g )+k(k + q)+m 2 1 
+ 2£ k (£ k - go ) 2 - #£ +q e^-iO + i 

1 £ k+ q(£ k+q + go)+k(k + q) + m 2 1 
2£ k+q (-Bk+ q + go) 2 - El e^+'j-") + 1 

1 £ k+q (£ k+q -go) + k(k + q) + m 2 1 
+ 2£ k+q (£ k+q - g ) 2 - El e«iWbO + i 

= ^[f(qo) + f(-q )]N F (E k ) (9) 
£ k (£ k -go) + k(k + q)+m 2 

(,-fck-go) -^ k+q 
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N F (E k ) - e0[Ek+ll) ~ 1 + ef}(Ek -n) + I 

For vacuum part where we obtained in last section, it describes the intrinsic graphene, in which the conduction 
band is empty while the valence band is fully occupied at zero temperature. When m — > graphene varies from a 
insulator to the zero- gap semiconductor system. When we take into account the finite density effect, the Fermi energy 
could lie either in valence band (/i < 0) or in conduction band(/i > 0) 



^matter -8 j ^L-[f (q ) + f (- q )]S (E k - fl) 



\Re ** f 1 ^(^-y^;/ - 1) (10) 



The retarded polarization function of the free fermions 



^ = \Re r ^ f ^ 4 ^-y k+ 2 g ° 2 - g2 %o 2 - 2^ k - - 2^) (11) 



7T 7o #k 7-1 \/T 



And in the above equations, we have defined 

Ref(qo) = [f(q )+f(-qo)]/2,ReS(f(q )) = [S(f(q ))-5(f(-q ))}/2 

If we set m = , we calculate explicitly and find they just coincide with the Ref.[4][5], and its finite temperature 
counterparts were discussed in Ref.[^||3| while we emphasize that within our model we do not need the overlapping 
factor used in their work and the results come out more naturally. 

SmC* tier = -Re f F dk f dx A ^' ~ 4k Qo±Qo ~ <f _ 2kqQ _ q 2 _ 2kqx ^ (12) 



7T Jo J -i yl - x 2 



^[ (g °: 2 2fc)2 „ 2 " 9 V {9( g - 9 o)e (fc -i±^) 
he(*-«)[0(^±5-*)-e(*^«-*)]} 

1 f** ^ r (g + 2fc) 2 -g 2 llr ^ g-gp 

— / d M 5 — 2 — ] 2 e(g-q )e(fc — — ) (13) 

n Jo q ~ go z 

fitter = ©(go - g)H+ + 6(g - go)n+ (14) 



Smn+ = -±=={{2k F - g0 ) Vg 2 - Wf - 9o) 2 + q 2 arcsin( 2 ^ g0 )6(g - \q - 2k F \) 

27r V ?o - g ^ 

2 

+ 7 ^-[Q(2k F -q -q)+Q(2k F -q + q)}} (15) 



o- tt+ 9(2fc F + g -g) r— — ■ -2 j 2 y/ (2fc F + g ) 2 - g 2 + 2fc F + g 

bmlTJ = {(2fc F + go)V ( 2fc f + go) 2 - g - g In 

27T V /g 2 - ^ g 
Um \ fjtyi V2 2 2, \/( 2k F ~ go) 2 ~ g 2 + 2fc F - ?0 Wol s-, 

- [(2fc F - qo)\J{2k F - go) 2 - g 2 - g In-* J6(2fc F - g - g)} 

(16) 
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The real part can be obtained using Kramers-Kronig relation 



Stolid = - r dq ^ mn f tteAq0) (17) 
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For the massive case(gapped graphene), the explicit calculation of the polarization function is straightforwardjl^ \h\ ■ 
so we just discuss an important static case qq = 



r* (q o) - — f kF ** T ix 1 f e ~ ^ 



k 



(\/^f + m2 — v m 2 ) 

2 f kF dkk 



irq Jo \Jk 2 + m : 



= (vy-4fc 2 - 4m 2 / vV - 4fc 2 )8(| - fc) 



7T 
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/ de(\Jq 2 + 4m 2 - 4e 2 



-Am 2 /y/q 2 + Am 2 - 4e 2 )9(V ' q 2 + 4m 2 - 2e) (18) 



2 /" M 

— / de(vV + 4m 2 - 4e 2 - 4m 2 /y/q 2 + 4m 2 - 4e 2 )6(vV + 4m 2 - 2e) 

1 .. g 2 - 4m 2 g ,_, , 

— { — mg H arctan Jfc)(2fc,F — gj 

7rg 2 2m 



n q 2 ~ 4m 2 J q 2 - 4k 2 F q 2 — Am 2 q , 

■ MV<T — 4fcp — mo arctan — — H arctan Bio — 2fcp) \ 

J v F 2 2[i 2 2to j v ;j 



IV. STATIC SCREENING 



(19) 



The static screening properties of the massive Dirac fermions in the RPA are controlled by the static dielectric 
function £rpa(o,0). 

27re 2 

e RPA (q,0) = 1 [n wac (q,0)+n matter (q,0)] (20) 

Kq 



Without the matter contribution, the polarized charge distribution has been discussed in Ref. 
account the matter part, things will be quite different 
For q < 2k F 



12j ,if we take into 



II(q,0) = -^ (21) 

7T 



Forg > 2k F 



2 M 1 ^, „, >, H, 777 q 2 -4m 2 Jq 2 - Ak% 



n(q,0) = — £ + — 9(g- 2k F )[fiJq 2 -4fc 2 - ^— arc tan^- ^] (22) 

it 15 v ^ 2 2/x 



Similar results has been obtained in Ref. 15]. From above we can deduce that the density of state at the Fermi surface 
is given by D{kp) = 2/i/tt, In the two limiting cases, the total static polarizability becomes a constant as in normal 2D 
electron liquid systems and 2D massless Dirac fermion systems. For 2D electron liquid, the 2D Thomas-Fermi wave 
vector is given by qpp = me 2 / k. Consider a external charge density embedded in the homogeneous fermions.Even 
though V e i(x) = C6(x) is short-ranged, the density response is not short-ranged which can be calculated by linear 
response theory Q. 

f2. 



H(X) = C / (2^ n(q '° )elq> 



(23) 



It is interesting to note that as the normal 2D electron liquid, the 2D Dirac fermions also has an oscillatory term due 
to the non-analyticity at q — hp. Next, we take into account finite temperature dependant screening, the polarization 
function at T ^ is 

n£ tter (q, T) = --{pL-m+\ HI + e-^-^} + ~ ln[l + e ~^+™)]} 
+ ^ r /2 de(lJt 2 0/2 -t 2 -2m 2 /Jel /2 -e 2 )[ ' ' ' 



nq 



9 /2 c ""'V «/2 n l + e^-») ^ l + e^+^r 

(24) 



At low temperature (T <C Tp),we have the asymptotic from of the polarizability from the above equation 

n(q,T)w-^p for e q / 2 < (J,, (25) 

In particular, for q = 2kp, m —> 0, we have 

n^ tter (q,T) = — — {ft + \ HI + e-^] + j ln[l + e"^]} 

4 f e "/ 2 I 1 1 



Trg J V 9/ 2 1 1 + ePk-ti 1 + e^+^ 1 
2 2k F f 1 

7T ' 7T 



M (T) + — L / dx^l - 2: 2 (1 - eM kFX -ri + e 2 ^*"^ + ...) 



ii 



And the total polarization function is 

II(q,T) « -1 M (D - 7X(1 - ^)C(|)T 3 / 2 (27) 
7T V ?rrc.F 2 2 

While for a general m, analytic expressions can not be obtained, however, it is quite easy to obtain numerical results 
from our semi-analytical results. 
At high temperature (T » Tp) 

n(q,T)«- M (T) + 4^-m 2 ] (28) 

For the chemical potential /Lt(T), we can get its explicit form the conservation of the total electron density. For the 
massless Dirac fermions, the asymptotic expression has been given in Ref. \v\ . We note that for the massive case the 
expressions is the same. 
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V. 



CONCLUSION 



In summary, we have presented a finite temperature field model for the 2D massive Dirac fermions and calculated 
the polarization functions for massive Dirac fermions. Finite temperature and finite density were taken into account 
to analyze the physics of gapped graphene in a more general case and connect various other group's important 



important quantities, for example conductivity or specific heat, can be obtain from our results. 

The polarization function at finite temperature can be also used to calculate the thermodynamic properties of 
massive Dirac fermion. 
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These results may be useful to study finite-temperature screening within the RPA. Other 



